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The cosmological peculiar velocity field (deviations from the pure Hubble flow) of matter carries significant 
information on dark energy, dark matter and the underlying theory of gravity on large scales. Peculiar motions 
of galaxies introduce systematic deviations between the observed galaxy redshifts z and the corresponding cos- 
mological redshifts z m ■ A novel method for estimating the angular power spectrum of the peculiar velocity field 
based on observations of galaxy redshifts and apparent magnitudes m (or equivalently fluxes) is presented. This 
method exploits the fact that a mean relation between z ras and m of galaxies can be derived from all galaxies 
in a redshift-magnitude survey. Given a galaxy magnitude, it is shown that the z cos (m) relation yields its cos- 
mological redshift with a lcr error of o\ ~ 0.3 for a survey like Euclid (~ 10 9 galaxies at z g 2), and can be 
used to constrain the angular power spectrum of z - z„(m) with a high signal-to-noise ratio. At large angular 
separations corresponding to / < 15, we obtain significant constraints on the power spectrum of the peculiar 
velocity field. At 15 < / < 60, magnitude shifts in the z m (m) relation caused by gravitational lensing magnifica- 
tion dominate, allowing us to probe the line-of-sight integral of the gravitational potential. Effects related to the 
environmental dependence in the luminosity function can easily be computed and their contamination removed 
from the estimated power spectra. The amplitude of the combined velocity and lensing power spectra at z ~ 1 
can be measured with < 5% accuracy. 

PACS numbers: 98.80.-k, 98.80.Es 



I. INTRODUCTION 

Dark energy and the theory of gravitation dictate the evolu- 
tion of large-scale structure in the Universe. The physical con- 
ditions allowing for the formation of galaxies, ultimately lead 
to a bias between the distribution of galaxies and the under- 
lying mass density. However, large-scale motions of galaxies 
are most certainly locked to the peculiar velocity field asso- 
ciated with the gravitational tug of the total underlying mass 
fluctuations. This assumes that gravity is the only relevant 
large-scale force and neglects the contribution from decay- 
ing linear modes. Despite the bias between the distribution of 
galaxies and the underlying matter field, the clustering proper- 
ties of galaxies have been the main tool for testing cosmologi- 
cal models. Currently planned galaxy surveys will allow us to 
quantify the clustering of galaxies on hundreds of comoving 
Mpcs and to even measure coherent distortions of galaxy im- 
ages which arise from gravitational lensing by the foreground 
matter. 

On the other hand, the peculiar motions of galaxies have 
traditionally been less successful as a cosmological tool, and 
there are several reasons for that. Neglecting other potentially 
important effects, peculiar velocities are approximately equal 
to the redshifts less the corresponding Hubble expansion re- 
cession velocities. The latter require direct measurements of 
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galaxy distances which are available for only a small fraction 
of galaxies. Presently, the number of galaxies with measured 
distances is several orders of magnitude below that of galaxies 
in redshift surveys used for clustering studies. Furthermore, 
although the underlying peculiar velocity is an honest tracer 
of the general matter flow, the inference of peculiar velocities 
from observational data is plagued with observational biases 
[Q]]. Traditional peculiar velocity catalogs are expected to im- 
prove within the next few years, but it remains questionable 
how well observational biases will be controlled, especially 
at large distances. An alternative probe of the peculiar veloc- 
ity field may be astrometric measurements of galaxies by the 
Gaia space mission [2, 3]. This probe is essentially free of 
the classic biases contaminating traditional peculiar velocity 
measurements, but it is also limited to nearby galaxies within 
~ 100/r 1 Mpc. 

Here we describe a method for deriving strong constraints 
on the power spectrum of the galaxies' peculiar velocity field 
independent of direct distance measurements and any bias- 
ing relation between galaxies and mass. The method is an 
extension of the approaches we have recently proposed |0- 
@], and it relies on using the observed fluxes of galaxies as 
a proxy for their cosmological distance [7]. Although this 
most basic "distance indicator" is very noisy, the large number 
of galaxies available in future surveys will allow one to beat 
down this noise to a sufficiently low level. Planned galaxy 
redshift surveys such as Euclid 0, [gl] will probe the structure 
of the Universe over thousands of (comoving) Mpcs, com- 
prising < 10 s — 10 9 galaxies at z ~ 1 and beyond. These 
observations will provide redshifts and fluxes, z and /, re- 
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spectively. The observed redshifts deviate from the cosmo- 
logical redshifts which would be observed in a purely homo- 
geneous universe. The large number of galaxies and the large 
sky coverage of these surveys can be used to derive a mean 
global relation between the mean redshift of a galaxy and its 
apparent magnitude m = -2.5 log/ + const. We interpret 
this relation as yielding the cosmological redshift z cos ('«) for 
a given apparent magnitude m. Angular power spectra of the 
difference zi - Z m {m.i) between the observed redshift z, of a 
galaxy and its expected cosmological redshift z cas should con- 
tain valuable information, mainly on the peculiar velocity field 
which is the main cosmological source for z, - z cm ( m d- In this 
work, we will show that the velocity power spectrum on large 
scales of a few 100 Mpcs could be constrained with signifi- 
cant signal-to-noise ratio (S IN) at the effective depth of the 
survey. Another contribution to z; - z cos (m,) results from the 
time evolution of the gravitational potential along the photon 
path, but is significantly smaller than that induced by pecu- 
liar velocities as we will show below. There are two addi- 
tional, indirect effects which modify the relation z cm ( m ) along 
a given line of sight. The first effect is related to the envi- 
ronmental dependence between galaxy luminosities and the 
large-scale structure in which they reside. Since this depen- 
dence is closely connected to the underlying density field, it 
can be self-consistently removed in our analysis. The second 
effect is caused by gravitational lensing magnification which 
changes the apparent magnitudes of galaxies in a given direc- 
tion. This latter contribution could actually be very rewarding 
since gravitational lensing provides a direct probe of the un- 
derlying mass distribution. Considering the analysis presented 
below, we will therefore treat it as part of the sought signal. 

The paper is structured as follows: We begin with a de- 
tailed description of the method and its application to galaxy 
redshift surveys in Sec. [TT] In Sec. [TTU we consider predictions 
for the standard ACDM model and discuss the method's via- 
bility as well as its expected performance. Finally, we present 
our conclusions in Sec. [IV] For clarity, some of the technical 
material is given separately in an appendix. In the follow- 
ing, we adopt the standard notation. The matter density and 
the cosmological constant in units of the critical density are 
denoted by Q and A, respectively. The scale factor a is nor- 
malized to unity at the present time (f = fo), and the Hubble 
function is defined as H = a/a. Further, r — c J'° df /a(t') will 
be the comoving distance to an object and z its corresponding 
redshift, assuming a homogeneous and isotropic cosmologi- 
cal background. Throughout the paper, the subscript "0" will 
refer to quantities given at t — to, and a dot symbol denotes 
partial derivatives with respect to time t, i.e. A = dA/dt. 



II. METHODOLOGY 

In an inhomogeneous universe, the observed redshift z of a 
given object differs from its cosmological redshift z cos (defined 
for the unperturbed background). The relative difference be- 
tween these redshifts, © = (z - z m )/(l + Z ras ) can be expressed 
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where r is a unit vector along the line-of-sight to the object. 
Here the radial peculiar velocity V and the usual metric po- 
tential O are assumed as relative to their present-day values 
at r = (/ = to). The first and second terms on the right- 
hand side of Eq. (Q~|i are the Doppler and gravitational shifts, 
respectively, while the third term describes the energy change 
of light as it passes through a time-varying gravitational po- 
tential. Note that this third term is equivalent to the late-time 
integrated Sachs-Wolfe effect experienced by photons of the 
cosmic microwave background (CMB). In what follows, we 
will denote the three terms as V , 0° and 0°, respectively. 
Also, we will consider angular power spectra (equivalent to 
angular correlations) of on large scales where the corre- 
sponding signal is significant only relative to the expected er- 
ror. Throughout this paper, we therefore rely on linear theory 
in a ACDM model where perturbations on all scales grow at 
the same rate. Introducing D(t) as the growth rate of the un- 
derlying mass density contrast 5 = p/p-1, linear theory yields 
the well-known relations 
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where £>( f o) = 1- The second relation is obtained from the 
first using Poisson's equation, i.e. V~0 = 3H l Q.od/2a. For 
the special case Qo = 1, we have D — a and fluctuations in the 
gravitational potential remain constant with time. 



A. Cosmological redshift versus apparent magnitude 

We aim to derive the field sampled at the positions of all 
galaxies in a flux limited redshift survey covering a significant 
region of the sky together with a large number of galaxies. As 
an example, we consider the planned Euclid redshift survey 
il|9l]. In order to obtain an estimate of 0, for each galaxy, we 
shall use apparent galaxy magnitudes (or equivalently fluxes) 
as a proxy to the cosmological redshift or distance. Although 
this most trivial distance indicator is very noisy, we will see 
that the high number of available galaxies allows one to beat 
down its scatter. 

The mean cosmological redshift Z ros (»i) corresponding to a 
given apparent magnitude m is 



z m ( m ) = z« 
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and the root mean square (rms) scatter around this relation is 
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where z\ and zi are the limiting redshifts in the survey, and 
P(z cos \m) denotes the probability that a galaxy with measured 
apparent magnitude m has a cosmological redshift z cos . The 
z cus ('«) relation can be linked to certain characteristics of the 
galaxy survey. The luminosity of a galaxy at z cos is L = A-nfd 2 L 
where di(z m ) is the luminosity distance to the galaxy, and its 
absolute magnitude is defined as M — -2.5 log L + const 
m - 5 log a!/, . Let <t>(M) be the underlying luminosity func- 
tion such that OdM is the number density of galaxies within 
the magnitude interval [M, M + dM]. Generally, the function 
O(M) depends on cosmic time, but in favor of a simplified de- 
scription, we will assume for the moment that it varies little 
throughout the depth of the survey considered. Note that it is 
trivial to include a time-dependent evolutionary term in <D(M). 
The probability P(z cm \m) satisfies 



P(z\m) oc P(m\z m )n(z), 
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where n(z m ) is the underlying mean number density of galax- 
ies at z m ■ In the absence of galaxy population evolution, we 
have 
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where d c is the comoving distance from the observer to z cos . 
The relation in Eq. (O can be easily derived using Bayes' 
theorem which yields 
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where P(m) is an observed quantity, P(z m ) is proportional to 
«(z cos ), and P(m\z cm ) is related to <t>(M) through 
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Here Mi = mi - 5 log djXz^,) is the absolute magnitude which 
corresponds to the limiting apparent magnitude m/ of the 
galaxy survey. 

Note that the scatter of z cos about the mean relation is not 
Gaussian. However, the associated lcr error in the derived 
angular power spectra depends only on the rms quantity <j z . 
Moreover, the central limit theorem implies that the errors in 
the power spectra tend to be Gaussian. Considering the actual 
observations, the two quantities z m im) and o~ z (m) can be esti- 
mated from the full survey by dividing the data in magnitude 
bins without actually computing P{z m \m). In this case, <T z (m) 
will include a positive contribution from the cosmological de- 
viations, but this is overwhelmed by both the intrinsic scatter 
in the z cos (»z) relation and the uncertainties in the photomet- 
ric redshifts. The underlying assumption is that the sought 
cosmological deviations between observed and cosmological 
redshifts cancel out when all galaxies of the entire survey are 
used. Clearly, a global constant mode should persist in this 
procedure, but we shall neglect it in this paper, assuming that 
it does not affect modes in the power spectra on smaller scales. 
Once we obtain z tos (m), we compute 
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for all galaxies in the survey. Note that in the denominator of 
the above, we have used z, instead of z cm (m). This substitution 
is consistent at linear order and motivated by the fact that z, 
is actually a better estimate of the true z m than z cm {m) which 
additionally includes deviations from the actual value of z cos 
due to large random errors with an rms of cr z . 

Environmental dependences of the luminosity distribution 
on the large-scale density field may systematically shift all 0, 
for galaxies lying in the direction of a certain line of sight. 
However, we will show below that the resulting signal con- 
tamination is small and in any case, it can be removed from 
the correlations since information on the underlying density 
field at the relevant scales will be directly available from the 
observations (see Sec. Ill Cl >. 



1. Expectations for Euclid 

As a test case for a future survey, we consider Euclid which 
aims at measuring photometric and spectroscopic redshifts of 
galaxies with z m ~ 1 over 15,000 deg 2 . More details on the 
Euclid mission can be found in 0. Here we focus on the 
photometric redshift survey since we aim at good statistics 
rather than precise redshift measurements. Photometric red- 
shifts will be measured with an error of cr phot < 0.05(1 + z P hot) 
for ~ 30 galaxies per arcmin 2 , with z tos > 0.7 and magni- 
tudes in the broad R+I+Z band (550-920 nm) RIZ AB < 24.5. 
Estimates of the photometric redshifts will rely on three near- 
infrared (NIR) bands (Y, J, and H in the range 0.92-2.0 fim) 
for objects with Y AB < 24, J AB < 24, and H AB < 24. These 
will further be complemented by ground-based photometry in 
the visible bands derived from public data or through engaged 
collaborations. 

The expected number density of Euclid galaxies with mea- 
sured photometric redshift can be parametrized as follows 
& 
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where zo = Zmean/ 1-4 12 is the peak of the redshift distribution 
and Zmean the median. Here we assume that z mea n = 0.9. We 
have compared the expected distribution n(z) of Eq. ( fT2b to 
that measured from galaxies with H AB < 24 and RIZ AB < 
24.5 in the zCOSMOS catalog (UHll. Indeed, we have 
found that Eq. ( TT2b does provide a good fit to the data in the 
range 0.7 < z p hot < 2.0 which we will consider in our analysis. 

The observed z p hot - H relation of zCOSMOS galaxies 
has been used to derive z p hot(#) and o- Zfhol (H) in different H- 
magnitude bins. To match Euclid constraints, we have only 
considered galaxies with RIZ AB < 24.5, Hab ^ 24, and 
0.7 < Zphot < 2.0. The additional constraints Y AB < 24 and 
J AB < 24 do not significantly modify our results and have 
therefore not been enforced. Our results for the zCOSMOS 
data are shown in Fig. Q] The (blue) solid curve and the (red) 
dashed curve represent the expected z m (H) and cr z (H) for Eu- 
clid galaxies, respectively. Since errors on the measured red- 
shift, i.e. cr phot , dominate the scatter for H< 20, they will be 
added in quadrature to determine the effective scatter in the 
Z m — H relation. 
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B. Spherical harmonics decomposition for discrete noisy data 

As usual, angular power spectra are defined in terms of the 
spherical harmonics Yi„,(r). For an all-sky continuous field 
f(r), the decomposition is 



flm 



-I 
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In our case, however, we have to deal with partial sky coverage 
(around 30% of the sky for Euclid) and consider that the field 
is sampled at discrete points given by the galaxy positions. 
The limited sky coverage could formally be described by an 
appropriate masking of the sphere ifl^flql . However, the de- 
scription in terms of masks will unnecessarily complicate the 
notation and somewhat obscure the physical interpretation of 
the results. We shall therefore resort to a simplified descrip- 
tion and assume that we are provided with a survey covering 
Anf sk steradians of the sky. For each degree I, we assume that 
there are (21 + l)f sky independent modes, instead of 21 + 1 for 
full sky coverage. Considering modes with angular resolution 
much smaller than the extent of the survey, we then have 
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where P/ is the Legendre polynomial of degree I. Throughout 
the paper, the angular integration is carried out only over the 



observed part of the sky, and the number of terms in the sum of 
the second relation is (21+ l)/ sk . As is obvious, these relations 
should be understood to hold in the approximate sense [17]. 

Considering a function f(r) with limited sky coverage in 
the continuous limit, we define 
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= fl! 2 f 



daf(f)Yi m (f), 
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where the integration is again taken over the observed region 
only. The angular power spectrum C/ is defined as the variance 
of the fi m 's and given by 



(16) 



where the average is taken over many different realizations of 
a field with the same power spectrum as /. For brevity, we 
will use the symbol < ■ > without any subscript to refer to this 
kind of averaging. Up to cosmic variance, the values of C; 
obtained from the ensemble average < \fi,„\ 2 > should equal 
those computed by averaging over the (21 + l)f sk modes for 
each /. Thus, the angular power spectrum may be estimated as 
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In the following, we shall interchange between these different 
kinds of averaging whenever appropriate. The factor 1 / f*' 2 in 
Eq. ( TT3T > and the rules in Eq. ( TBI guarantee that 
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which can be shown by decomposing Pi into F/ m according 
to the second relation in Eq. ( fPfl i and assuming that C{9) is 
negligible for angular separations larger than the extent of the 
survey (see AppendixlAl. 

For a discrete sampling of / at the positions of N galaxies 
distributed over the observed part of the sky, we write ifTill 



N „ 

2] Ah)* J n(r)f(f)dQ. 
*hj f(r)d£l, 



(19) 



where n{r) is the projected number density of objects and h = 
N/(f sk An) is the corresponding mean number density over the 
observed part of the sky. Here we have assumed that / itself 
depends on the density contrast 8 — n/h— 1, so that the last step 
applies to linear order in the fluctuations. Again, the angular 
integration is carried out only over the observed part of the 
sky. In analogy to Eq. (15[ , we further define 



galaxy according to the cr, in the sum of Eq. d20l i. To min- 
imize the effects of shot noise, we weight each galaxy by a 
factor of Wi which is given by 



No-T 



wj = 



This particular weighting scheme yields 



N 



-2 ■ 



(26) 



(27) 



where we have assumed that that the underlying signal makes 
a negligible contribution to cr 2 .. The weighting does not af- 
fect the ensemble average of |//„,| 2 , and its net effect is that 
cr/ should be computed as given by Eq. ( f27l ). i.e. it merely 
reduces shot noise errors (for additional details, see Appendix 
IBl . In principle, one may use any weighting scheme. 
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In this case, one can show that the angular power spectrum 
takes the form (see Appendix iBb 



Cj =< l//„,| 2 > ~, 
11 



(21) 



where the second term represents the contribution of shot 
noise due to the discrete sampling of / and cr 2 f = /(^V) 2 IN. 



To estimate the values of C/, we therefore use 



C, =< \f„„\ 2 > - 



(22) 



If fj — S (r,) + e, where S is an underlying cosmological signal 
and e, is an uncorrected random error, this simply reflects the 
fact that Ci =< \S i m \ 2 >. The expected error in this estimate 
of Ci is given by (again, see Appendix IBl 
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which includes contributions from both shot noise S sn and 
cosmic variance Z cv , 



The variance of the scatter <rf, 
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depends on both magnitude and redshift. Note that the factor 
of (1 + zd arises from the definition of in Eq. (Q]). Consid- 
ering the application to real data, it is prudent to weight each 



C. Environmental dependences in the luminosity functions 
and magnification by gravitational lensing 

So far, we have assumed that the systematic shifts in the 
Z m ( m ) relation for galaxies in a given direction are solely due 
to the terms appearing on the right-hand side of Eq. ([1]). How- 
ever, additional shifts may arise from changes in the mean 
magnitudes due to large-scale density fluctuations in a given 
direction, i.e. environmental dependences and evolution in the 
luminosity function, and the magnification effect caused by 
gravitational lensing. In what follows, we shall denote their 
contribution as env and © lens , respectively. Both effects re- 
sult in a magnitude shift which translates into additional cor- 
related deviations of the estimated z an from its mean relation. 
To model the contribution of these effects in the correlation 
function of 0, we need to translate a magnitude shift Am into a 
corresponding shift Az. This can directly be read off the (blue) 
solid curve z m (m) shown in Fig. [1] in which a shift Am = 0.2 
leads to mean shift Az Lm ~ 0.017 IT19I1 . 

To quantify the impact of environmental dependences, we 
assume that the systematic shift in galaxy magnitude depends 
on the density contrast 5 in regions where they reside, irre- 
spective of the smoothing scale l20l l2lll . Here we adopt the 
linear relation Am = 0.26 which is observationally inferred 
from the visible band 12211 . The actual dependence is a func- 
tion of the photometric band, and it is expected to be much 
weaker in Euclid's H-band as indicated by the weak environ- 
mental dependence of the Schechter parameter M* fitted to 
NIR luminosities 12311 . The remaining Schechter parameter, 
usually dubbed a, exhibits a stronger dependence, but since 
this parameter fixes the shape at the faint end, it will have 
only a small impact in deep surveys like Euclid. Nonetheless, 
it will be possible to remove most of the contamination caused 
by this effect using the observed distribution of galaxies in the 
Euclid survey. 

In the weak-field limit, the magnification induced by grav- 
itational lensing is proportional to 1 + 2k, where k is the ef- 
fective convergence field, i.e. an integral over the (weighted) 
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density contrast along the line of sight 12411 . For a fiat ACDM 
model and a fixed source redshift corresponding to a comov- 
ing source distance r, one finds 



k(9, r) = 



3// 2 Qq 

2c 2 
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Jo 



r'(r - r') 

dr'— -6(r'e,r'), 

a(r')r 



(28) 



where the two-dimensional angular vector 9 is perpendicular 
to the line of sight. Therefore, it is straightforward to model 
the magnification and its effect in our analysis. The magnifi- 
cation field contains valuable information since it probes the 
growth of the angular derivatives of the gravitational poten- 
tial fe. g. I24L I25I1 . It can also be used to constrain the gravita- 
tional slip which arises in certain modifications of the general 
theory of relativity [e.g. H^, [27l]. Given a theory of gravity, 
much of the contribution to the power spectra from gravita- 
tional lensing can, in principle, be removed using the under- 
lying large-scale density field which is inferred from the fore- 
ground galaxy distribution. Considering the following analy- 
sis, however, we will treat the effects of lensing magnification 
as part of the signal. As we will see below, its contribution is 
not negligible, and can be constrained together with the power 
spectrum of the velocity field. 



III. THEORETICAL ANGULAR POWER SPECTRA 

In the following, we will present predictions for the angular 
power spectra corresponding to the various terms in Eq. (Q]). 
Since the observed galaxies cover a redshift range z\ < z < Zz, 
we do not consider the angular power spectra of a quantity / 
defined at a specific redshift, but instead we use 



f[rr,t(r)] p(r)dr, 



(29) 



where r\ and rj_ are the comoving distances at the survey's 
limiting redshifts, z\ and zi, respectively, and p(r)dr is the 
probability of observing a galaxy within the interval [r, r+dr]. 
Note that the function / represents any of the terms on the 
right-hand side of Eq. ([T), i.e. V , 0°, and 0°. Beginning 
with 0°, we write 



0* /m = J" dD^YUf) 
= i J" dClY lm (¥) J 



(30) 



dClYUr) Wdr)® (rr)dr, 



where we have used the linear relation 0(r, t) = (D/a)(!>o(r, to) 
and defined W$ = Dp(r)/a, with D(t) and a(t) evaluated at 
t — t(r). Expanding <l>o(r) in Fourier space, 

®° (r) = J dh®^, (31) 



and using 



An^i'MkrWimYUk), 
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where ji is the usual first-kind spherical Bessel function of 
degree I, we get 



<2> lm 



2^? J 2 drW * f d * k(S) kj' (k ^ Y t»'&- ( 33 > 



Therefore, we finally arrive at 



Cf =< |0°/ m | 2 > 



J dk^P^k) f drWfjiikr) 
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where we have used < fl^ib' >= (2n) 3 6 D (k - k')P®(k). Sim- 
ilarly, using the linear relation in Eq. ©, we obtain 



C 



j dkk 2 P^(k) J " drWy |y 



kj 
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where W v = 2aDp(r)/3Q ( )H ( 2 1 . As for © 6 , the last term ap- 
pearing in Eq. (Q~|), we will assume that the signal is mostly 
caused by the large-scale structure between the observed high 
redshift galaxy sample and the observer. Therefore, 0° is ap- 
proximately the same for all galaxies along a common line-of- 
sight such that 0° * 0°. The angular power spectrum then 
reads 



Cf =< |0\,| 2 > 



J dkk 2 Pq,{k) j dr'W^ji{kr') 



(36) 



where Wj,(t)/a = (d/dt)[D/a]. The integration over r' is taken 
from r — out to a distance beyond which <S> becomes nearly 
constant with time. For simplicity, we will assume that most 
of the contribution to the integral comes from the inner edge 
of the considered galaxy sample, and adopt z = 1 in Eq. (|36] l 
for all galaxies in the survey. For Euclid, this seems to be the 
case since the mean redshift is expected to be z ~ 1 . Note that 
the same assumptions are used when calculating the angular 
power spectrum C lens . 



A. Predictions of the ACDM scenario 

Having derived the relevant expressions above, we are now 
ready to make predictions for the framework of ACDM. Here 
we adopt a spatially flat model with best-fit parameters based 
on the CMB anisotropies measured by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 028I1 . In this case, the 
total mass density parameter is Q,„ = 0.266, the baryonic den- 
sity parameter Qb = 0.0449, the Hubble constant h = 0.71 
in units of 100 kms^'Mpc -1 , the scalar spectral index n s = 
0.963, and erg = 0.80 for the rms of linear density fluctuations 
within spheres of 8/j~' Mpc. We work with a parametric form 
of the power spectrum taken from Ref. lf29ll (see Eqs. 29-31 
in their paper). For the calculation of angular power spectra, 
we use the expressions from Sec. [Ill] together with p{r) cor- 
responding to the redshift distribution of galaxies appropriate 
for Euclid which is given by Eq. ( fT2] i. 
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FIG. 2. Angular power spectra in the ACDM model as explained in the text: The shaded area represents the cosmic variance uncertainty E cv 
on C tot , and E sn is the lcr shot-noise error for the Euclid survey. The desired signal is contaminated by C cnv which results from environmental 
dependences in the luminosity function and can be removed from the actual Euclid data with high precision represented by E cnv . Only 
contributions from V (peculiar velocity), 0* (gravitational shift), and lells (lensing magnification) are included in C tot . At / ~ 60 (not 
shown), S sn and C tot become comparable. 



The angular power spectra and associated errors are plotted 
in Fig. [2] The shot-noise error S sn is computed using Eq. d27l i 
and assuming crj = 0.17 as inferred from the (red) dashed 
curve in Fig. Q] For completeness, we also show the power 
spectrum C env which results from environmental dependences 
in the luminosity function. As explained in Sec. Ill CI the 
contamination arising from this effect can be removed given 
the observed distribution of galaxies in the survey. The accu- 
racy to which this can be achieved for Euclid is represented 
by 2 e nv, the If error within which C env can be estimated ex- 
plicitly from the data (see Appendix [B] for details). As is 
clear from C lens in the figure, magnification caused by grav- 
itational lensing introduces significant angular correlations in 

(Zi-Z c J/(l+Zi). 

The quantity C tot is the power spectrum of the sum of the 
three signals © lens , ® v , and ®°, corresponding to the lensing 
magnification, the Doppler shift, and the gravitational shift, 
respectively. Note that the calculation of C tot does include 
covariance between the three individual signals. The contri- 
bution from ®° is negligible as is indicated by C°, and we do 
not include it in C tot . At / £ 10, C v is the dominant contri- 
bution, but the lensing term C lens takes over at I > 10-15, 
roughly until / ~ 60 where it drops below the shot-noise level. 
The shaded area represents the lcr error of C tot due to cos- 
mic variance. For practical purposes, it is therefore possible 
to provide measurements of C v and C lens by an appropriate 
fitting procedure of the two corresponding curves to the mea- 
sured C tot . Unfortunately, at low I, cosmic variance is so large 
that accurate constraints on C° do not seem possible. 



B. Signal-to-noise and expected error on model normalization 

Suppose that observations yield an estimate C° bs for the to- 
tal power spectrum Cf l . Let Cf° = and Cf' be the expected 
total power spectra for the null hypothesis with no correlations 
(Ho) and of the ACDM model (Hi), respectively. The ratio of 
probabilities for H\ and Hq is given by 



-2\nP(H )/P(Hi) = 



S 



(fobs — C H °\ (fobs _ f H i\ \ 



(37) 



If H\ holds, it follows that C, 1 is equal to the C ; tot shown in 
Fig. H and C° bs - Cf is a random variable with variance 2 
given by Eq. (f23t . Therefore, the signal-to-noise ratio (S IN) 
for rejecting Ho is 



(||J = -2\nP(H )IP(H x ) = 

r i 



(2/+D/ sty C ; 2 



nCi 



-21n 



1 + 



nCi 



(38) 



where Ci = C, 



^1 = o, - = Cf 1 . Substituting the relevant quantities 
into the above, we obtain (S /N)\o = 101, where the above 
sum rapidly converges by / = 30. On the other hand, if Hq is 
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i 1 1 1 1 1 1 r 




I i i i i i i i i I 

0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 

FIG. 3. Confidence levels for measuring the amplitudes of the velocity and lensing power spectra by fitting aJC v + c^ ms C ,sns to random 
realizations of C tot which include both shot noise and cosmic variance: The fitting procedure ignores physical cross-correlations between the 
used quantities. Contours are shown for A^ 2 = 2.3, 6.17, and 9.2, corresponding to confidence levels of 68%, 95.4% and 99%, respectively. 



true, then the S /N for rejecting Hi is given by 



two-parameter model 



s 



-21nP(// 1 )/P(//„) 



(2/+D/ sky C ; 2 



Kaj/n + Ci) 2 



1 + 



nCi 



• f ' 



+ 21n 



1 + 



nCi 



(39) 



which quickly converges at I ~ 50 and yields (S /N)\\ = 14.7. 
The fact that (S /N)\\ <k (S /N)\o is a result of cosmic vari- 
ance which is zero for the Ho hypothesis, but significant in the 
ACDM scenario, i.e. £ 2 V = 2C / 2 /(2Z + l)/ sty . 

In addition to S /N considerations, we can assess how well a 
measurement of C tot would constrain the normalization of the 
ACDM power spectrum in terms of <x g . To this end, we write 
the model's total power spectrum as C tot m = (cr 8 /0.8) 2 C tot , 
where C tot as illustrated in Fig. [2]is obtained for erg = 0.8. The 
expected Icr error on erg is then (-d 2 \nP(Hi) I dcrfy~ l , where 
P{H\) is now expressed as 



lnTO)^-^ 



(c° bs - c; o, ' m ) 2 



2I 2 



(40) 



Cf t,m Using a normal 



and £/ is given by Eq. d23l l with C/ - 
probability distribution of the form given in Eq. ( 1401 with 
fixed cr 8 = 0.8, we have generated 1000 random realizations 
of C obs . For each of these realizations, we have maximized 
P{H\) given by Eq. (l40l with respect to erg. As a result, we 
find that the true value erg = 0.8 is recovered within a relative 
lcr error of less that 4%, without any statistically significant 
bias. 

We have also inspected the possibility of constraining the 
velocity and lensing signal amplitudes with the help of the 



C tot ' m = a Y C 



V 2 

+ a 



lens 



c 



lens 



(41) 



This model for C tot m neglects the contribution of 0° to C obs 
as well as any covariance between the remaining signals, V 
and lensing magnification lens . Using the above expres- 
sion, we have repeated the procedure described above for 
constraining erg with the 1000 random realizations of C obs . 
The result is presented in Fig. [3] which shows contours of 
A^ 2 = -21nP + 21nmax(P) as a function of ay and ai en s- 
The contours have been computed for one of the random re- 
alizations, but the values of the best-fit parameters (giving the 
lowest A^ 2 ) obtained for this particular realization have been 
shifted to their underlying value, i.e. unity for both. This is 
reasonable since the average best-fit parameters of 1000 real- 
izations are essentially unbiased. As seen from the figure, ai ens 
is constrained with better accuracy than the velocity ampli- 
tude ay. This is consistent with Fig. Q]which shows that C lens 
dominates the total signal over a large range of / while C v is 
significant only for low / where cosmic variance becomes in- 
creasingly important. Still, both parameters are constrainted 
with good accuracy. 



IV. DISCUSSION 

In this paper, we have presented a novel method for de- 
riving direct constraints on the peculiar velocity and gravita- 
tional potential power spectra from currently planned galaxy 
redshift surveys. The large number of galaxies with photo- 
metric redshifts in these surveys allows one to exploit appar- 
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ent galaxy magnitudes as a proxy for their cosmological red- 
shifts since it beats down the large scatter in the z On) relation 

o COS v ' 

and the uncertainty in the photometric redshifts. The method 
aims at directly constraining power spectra of the underlying 
fluctuation fields independent of the way galaxies trace mass. 
Other methods for extracting cosmological information from 
redshift surveys rely on accurate measurements of the galaxy 
power spectrum in redshift space (since galaxy distances re- 
main unknown). The power spectrum and other statistical 
measures based on the distribution of galaxies have been suc- 
cessful at probing the nature of dark matter and placing im- 
portant constraints on neutrino masses |[30l - [32tl . Having said 
that, however, they depend on a very accurate knowledge of 
the relation between galaxies and the full underlying matter 
distribution. The method we have considered here is less pre- 
cise, but it is completely independent of the galaxy formation 
process and offers a much more sensitive assessment of the 
underlying physical mechanism driving cosmic acceleration 
and structure formation. This approach is particularly worth- 
while if such constraints on the velocity field and the gravita- 
tional potential are contrasted with local constraints obtained 
from data at low redshifts. For example, peculiar motions of 
galaxies within a distance of ~ 100/i~' Mpc can be measured 
usin g tig ht relations between intrinsic observables of galaxies 
le.g |33L l34ll . and also using astrometric observations of the 
Gaia space mission which is currently scheduled for launch 
in 2013. These peculiar motions of galaxies have been useful 
for constraining cosmological parameters l35ll as well as the 
amplitude of the velocity field in the nearby Universe JH [3(1 . 

Although we have presented predictions for the Euclid sur- 
vey, the science proposed in this paper will not have to wait 
for this space mission. In fact, several ground-based photo- 
metric surveys in the optical and near-infrared bands, which 
will constitute the backbone of Euclid's photometry, will pro- 
vide photometric redshift catalogs that can be used for our 
purposes well before the launch of the satellite. On a shorter 
timescale, the VLT Survey Telescope (VST) will be used to 
carry out the Kilo Degree Survey (KiDS), one of the ESO 
public surveys. It will cover 1,500 deg 2 to u = 24, g = 24.6, 
r = 24 A, and j = 23.1, and will probably contain ~ 10 8 galax- 
ies with measured photometric redshifts ll37tl . Also, the Dark 
Energy Survey (DES) will start its operations soon, and it will 
cover 5000 deg 2 of the Southern sky within 5 years, reaching 
magnitudes up to ~ 24 in SDSS griz filters, comparable to the 
limiting magnitudes of Euclid and with a redshift distribution 
dN/dz similar to that of Euclid galaxies. DES will measure 
photometric redshifts of ~ 3 x 10 s galaxies with cr photo ~ 0.12 
at z ~ 1 13811 . Furthermore, the first of the four planned 
Pan-STARRS telescopes l39ll has been operational since May 
2010. The planned 3n area of the sky will be considerably 
shallower, detecting galaxies below a limiting magnitude of 
~ 24 in the griz bands. A deeper survey involving the PS 1 and 
PS2 telescopes is currently being planned. The survey should 
cover ~ 7, 500deg 2 with limiting fluxes g = 24.7, r = 24.3, 
i = 24.1, and z = 23.6. Photometric redshifts will then be 
measured for ~ 4.5 x 10 8 galaxies with similar errors. Fi- 
nally, on the long run, the Large Synoptic Survey Telescope 
(LSST) is expected to start operations in 2020. Its main deep- 



wide-fast survey is expected to observe ~ 20, 000 deg 2 in the 
ugrizy bands. After about 10 years of operation, it will reach 
much deeper depths (down to a co-added magnitude r = 27), 
detecting about 3 x 10 9 galaxies l40ll . 
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Appendix A: Proof of the relations in Eq. fl8t 

For any function /(r), we have defined 



fim = 4/2 f dnf(r)Y lm (r) ■ 

/sky J 

Using the relations given in Eq. ( fT4b . we get 



(Al) 



C ' = (2/ + 1 l)/ sky ? <l/ ' m|2 >= (2/ + W 2 y 
x f dQdQ! < /(r)/(r') > £ Y lm {r)Y;, n {r') 

= —5— f dQ.dQ.'C(r ■ r')P,(r ■ r') 

= f dC1 ' f d£lC(r-r')P,(r-r') 



(A2) 



c/QC(cos 6>)P/(cos 8), 



where for the last step, we have used that the integral over 
dQ. is independent of r', assuming that the coherence angular 
length is smaller than the spatial extent of the survey. Sim- 
ilarly, the first relation in Eq. ( fTSl follows from substituting 
it into the second and exploiting the (near) orthogonality of 
Pi(cos0). 



Appendix B: Estimating angular power spectra for discrete data 

Just as in the main body of the paper, we will take the sym- 
bol < (■) > x = (1/(2/ + l)/ sty )2m(0 to denote averaging over 
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all indices m corresponding to the degree / of the spherical 
harmonic decomposition. Again, the symbol < • > without 
any subscript will refer to the previously introduced ensemble 
average, and we shall interchange between < • > m and < • > 
whenever appropriate. For brevity of notation, we also use 
the definition Y' m = Yi m {ri), thus omitting the subscript / since 
all calculations will refer to a given degree / of the spherical 
harmonics. Starting from the definition of fi m in Eq. ( l20t . we 
obtain 

^ ' m ' m 



< \flm\ 2 >= Z^T- J] < fifj > YirJ 



' sky 



n 



(Bl) 



sky ft j 



where the value of err is inferred from the z. Am) relation, i.e. 
cry = Yii °i IN- As our estimate for C/, we therefore take 



Ci =< \f lm \ 2 > 
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(B2) 



Another way of arriving at this result is to partition the ob- 
served sky into infinitesimally small cells of ang ular size £Q, 
with each cell containing at most one galaxy Illql . In this case, 
we write 



flm — _ 1 12 / j fa n a Y m , 



(B3) 



where n a , the number of galaxies in cell a, is either or 1 . 
Using that < n a >- hdQ. in this representation, we obtain the 
expression 

1 



< l//ml > 



" 2 / s 



aWft - > Y m Y m 



n 2 L 



n 2 L 



Y J <f 1 a ><nl>\Y a „\ 2 +C l 

a 

Y,<fl>n5Cl\Y<,;\ 2 + C l 



(B4) 



= -r+C u 
n 

where we have used the relation n 2 , = n a . We find this partic- 
ular approach a little more intuitive to compute the expected 
error of the estimate in Eq. (IB21 >. Let f and n r be possible 
random realizations of the data in another universe. The vari- 
ance of the error in C/ may then be written as an ensemble 
average over these realizations, i.e. 

s 2 = < (cr - c,) 2 > 



i 



« 2 / s 



1 



^jfafp n a n f} < Y m Ym 



> -Q 



h 4 f 2 



2 < f> r jw> r «'fh nr v 



(B5) 



^ 'm'm > m < ~ 1 m' 1 m' > ■ 



Note the remarkable fact that due to the uncorrected nature of 
the errors and because a + ft as well as a' + ft' , only second- 
order moments of « and / contribute to S 2 . This is important 
since it implies that the non-Gaussian nature of the scatter in 
the z ctls {m) relation does not affect the variance of errors in the 
estimated C/. Assuming that fn is a Gaussian random field, 
we further have 



(B6) 



= < fafpn>0 x fZfpt&ifp > 

+ < faf»'a> >< ffifptfrf, > 

+ <£/W 

Considering the above, averages with no equal indices will 
contribute terms proportional to Cf in the sum of Eq. ( IB5I ). 
Neglecting the contribution of shot noise, this would yield the 
usual cosmic variance expression 2Cf/(2l + l)f sky - For the 
case a - a' and /3 = /3' , the second term on the right-hand 
side of Eq. dB6b turns into <r^(h6Q.) 2 , where we have again 
used that < n 2 >=< n a >= hSQ.. The other combination 
which makes the same contribution is a — /?' and a' = ft. 
A bit of algebra shows that the contribution to the variance 
in Eq. JB5I ) of these two combinations is 2o"i/(2Z + l)« 2 / sty . 
For a - a' and ft + ft' , the second term on the right-hand 
side of Eq. dB6l ) simplifies to cry < fjjfp, > («5Q) 3 , and the 
total result of similar permutations with only two equal indices 
reads 4cr 2 C//(2Z + l)«/ sky . Finally, summing up all relevant 
terms leads to 



E 2 = 



(2/+D/ sk 



— + Q 

n 



(B7) 



Using the weights defined in Eq. ( 1261 ). it is easy to say that the 
same results are valid with cr y given by Eq. ( l27l i. 

In the main body of the paper, we have assumed that en- 
vironmental dependences of the galaxy luminosity function 
arise due to variations in the large-scale density contrast. The 
signal contamination associated with this effect has an angu- 
lar power spectrum which is proportional to that of the density 
contrast. In the following, we give an estimate for this power 
spectrum and the Icr error within which it can be constrained 
from the observed galaxy distribution. To begin with, we con- 
sider |E] 



Sim - 



1 



- All 

nf4 



(B8) 



For the fi m 's, we have assumed that n a is a discrete sampling 
of a uniform distribution. That was consistent with linear the- 
ory since we have assumed that / is a function of the density 
contrast. However, now we are interested in the autocorre- 
lation function of n a , and thus the treatment of <5/,„ is a little 
different. To obtain an estimate for the angular power spec- 
trum of the density contrast, we write 



<\S b 



(B9) 



= -+c°. 

n 
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Similar as before, we hence take Cf =< |<S/ m | 2 >„, -\jn as our 
estimate of the angular power spectrum of 5. The statistical 
error in this estimate can be derived using the same approach 
as above, and it is given by 



E-L = (- + Cf) . (BIO) 



One of the main interests in this paper is to remove the con- 
tamination C™ v to the signal Cf 1 using the actually observed 
distribution of galaxies. Note that in this case, cosmic vari- 
ance is irrelevant to the analysis since we are referring to 
a given data set rather than the expected deviation of the 
estimated power spectrum from some underlying theoretical 
model. Therefore, we may simply neglect the contribution of 
cosmic variance in Eq. ( IB 101 ). and the variance of the relevant 
error reduces to the expression 2/(2/ + l)n/ skv . 
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